Add primerprospector_analyzeprimers#11615
Conversation
timrozday-mgnify
left a comment
There was a problem hiding this comment.
Thanks @vagkaratzas ! This looks useful.
I have one major issue with this module relating to incorrect handling of output hit files when there are multiple primers in the input primer file (a common case). The tests should be amended with this case.
I think it would be useful to support gzipped fasta inputs, but I don't know what nf-core good practice says about that.
| -f "${fasta_arg}" \\ | ||
| ${primer_arg} | ||
|
|
||
| mv *_hits.txt "${prefix}_hits.txt" |
There was a problem hiding this comment.
For multiple primers in the input primer file this will cause an error or overwrite some results. Rather all the files should be emitted in the hits output channel. If prefix renaming is desired then it needs to be performed on each of the produced hits files.
There was a problem hiding this comment.
Good catch, updated!
Since the module allows a list of fasta files (and not just a single fasta one), it is a bit more complicated to do that in the code (but definitely not impossilbe). But for now I guess running |
timrozday-mgnify
left a comment
There was a problem hiding this comment.
Looks great, thank you @vagkaratzas!
Just two edits to the error messages since prefix is no longer used to rename output files.
I've approved so you can crack on with merging when you choose.
Co-authored-by: Tim Rozday <timrozday@ebi.ac.uk>
Co-authored-by: Tim Rozday <timrozday@ebi.ac.uk>
* bots initial attempt * removed unnecessary temp output dir, and renamed output files, added comment for bug workaround * flagged the bug in the comment * execfile on PATH for shim, instead of local user/ path * stricter version eval to pass GitHub runners warnings * updated and tested for multiple primers * Update modules/nf-core/primerprospector/analyzeprimers/main.nf Co-authored-by: Tim Rozday <timrozday@ebi.ac.uk> * Update modules/nf-core/primerprospector/analyzeprimers/main.nf Co-authored-by: Tim Rozday <timrozday@ebi.ac.uk> * prefix only used in stub --------- Co-authored-by: Tim Rozday <timrozday@ebi.ac.uk>
Workaround shim in main:
PR checklist
topic: versions- See version_topicslabelnf-core modules test <MODULE> --profile dockernf-core modules test <MODULE> --profile singularitynf-core modules test <MODULE> --profile condanf-core subworkflows test <SUBWORKFLOW> --profile dockernf-core subworkflows test <SUBWORKFLOW> --profile singularitynf-core subworkflows test <SUBWORKFLOW> --profile conda